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Variable Structure PID Control to Prevent Integrator Windup 
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Figure 1: PID control with saturation limits 

1 Introduction 

PID controllers are frequently used to control systems re- 
quiring zero steady-state error while maintaining require- 
ments for settling time and robustness (gain/phase mar- 
gins). PID controllers suffer significant loss of perfor- 
mance due to short-term integrator wind-up when used 
in systems with actuator saturation (see Figure 1). We 
examine several existing and proposed methods for the 
prevention of integrator wind-up in both continuous and 
discrete time implementations. We may write a continu- 
ous time PID control law as 

K(s) = ^1 = K P + sK d + Kr/s (1.1) 

where u n (t) is the nominal control command and e(t) = 
yref{t)~y{t) Is the error between a reference signal y re /(£) 
and output y(t) of the system being controlled. The re- 
spective state-space implementation is 

fj = e u n — Kpe + Koe + K/t] 

Control saturation occurs when u n lies outside of actuator 
limits, u n ^ [u m i n ,u m ax]* 

2 Background 

Although numerous methods have been proposed for the 
prevention of windup in controller integrators and slow 
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dynamics, very few’ textbooks discuss the problem, par- 
ticularly at the undergraduate level (e.g., [9]). The earliest 
treatment of anti-windup techniques that w^e are aware of 
was done by Fertik and Ross [8]. This technique fits into 
the larger class of anti-windup bumpless transfer (AWBT) 
control, example methods of w’hich are covered in [1], [2], 
[3] , [4], [6], [11], [12], [13], [14], [15], [18], [22], [23]. A 
general theoretical framework for the parameterization, 
synthesis and analysis of AWBT control is provided in 
[18], in which it is assumed all nonlinearities are exter- 
nal to the controller; i.e., the controller is required to be 
linear. Stability analysis for these methods is typically 
performed through describing functions, see, e.g., [2]. 

Other anti-wdndup methods include conditional integra- 
tion and/or integrator limiting, (e.g. [7], [10], [16]) which 
freezes or “clamps” the integrator value when certain con- 
ditions are not met, e.g., saturation, output not in “pro- 
portional band,” etc., the use of time varying gains to 
avoid saturation [17], or the setting of the integrator to a 
prescribed value during saturation, also called preloading 
[21]. These methods do not fall into the class of AW BT 
control, since the switching action on the integrator ren- 
ders the method nonlinear. 

3 VSPID control 

While many of the above methods are applicable to mul- 
tivariable systems, w'e shall confine our attention in this 
paper to the treatment of PID control with individual 
saturation limits for each PID channel. W T e shall contrast 
three methods of conditional integration (Cl), one method 
of “preloading,” a simple AWBT method, and a new vari- 
able structure PID (VSPID) controller. Our discussion 
makes use of the following definition. 

Definition 3.1 The saturation function 

sat(a, a m i n , a max ) = max(u m j n , min(a, a maa: )). 

It will be seen that the VSPID and AWBT methods yield 
similar behavior when the AWBT uses a high-gain feed- 
back of the control saturation error u n — u where we 
define 

u s = sat(u n , u m j n , u max ) (3-2) 






Tom porn turoi 


Methods that we examine here are: 

CI-I Integrator limiting; see [5], p. 278. Impose hard 
limits (saturation) on the integrator value q: 

{ 0 T) i [r)m,n,Vmax] and e X (T) - f}) > 0, 

— {jjmin d - r7 max )/2 

e otherwise 

The choice of design parameters q mi n,qmax is not always 
clear; for this study, we choose 

(jlmimVmax) — 

CI-II Freeze integrator input q at 0 when u n is in sat- 
uration: 

. _ f 0 u n ^ u s (see Equation (3.2)) 

^ ~ | e otherwise 



Figure 2: Closed loop temperature profiles of simulated 
furnace with continuous time anti-windup PID feedback 
laws 


CI-III Freeze q when u n is being driven into saturation; 
that is, 

. f 0 Un and e(un — u$) 0 

^ ~ \ e otherwise 

(Kf > 0 is assumed.) 

Preloading Manually reset integrator value q to an off- 
line predetermined value qd when u n is in saturation. For 
the purposes of simulation, we implement this technique 
by modifying the integrator input as 

. _ f -a x (q - q d ) u n ^ u 9 
^ 1 e otherwise 

where the parameter a > 0 controls the integrator decay 
rate when u n is in saturation. For this study, we select an 
integrator value of q d = 0. 

AWBT [9] p. 198. Include an integrator feedback term 
in the integrator involving the error between the nominal 
control u n and its limited value u s : 


VSPID Rather than freeze the integrator value as shown 
above, dynamically drive the integrator so that u n lies at 
the edge of the saturation region: 

r Un ^ Us and e(l 77 u) - 55 0< 

’j = 1 u = {Umin + U max )/2 (3-4) 

I e otherwise 

where a > 0 is a positive constant selected such that u n 
rapidly converges to the nearest extreme value of [u min > Umax 
Discrete time implementation of these control laws is 
straightforward. 


4 Simulation examples 


The anti-wdndup methods of the previous section were 
simulated in closed loop with a model of an electric fur- 
nace P(s) with state-space model 
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i;(t) (4.1) 


q — 6 R"a(^n ^s)* 

Notice that u n ^ u s implies that the error e and the 
additional feedback term must “fight” one another. This 
property is inherent due to the linear nature of the control 
law. 


where v(t) is an input voltage constrained between 0V 
and 10V, f(t) is the filament temperature and c(t) is the 
chamber temperature. The uncompensated settling time 
of the system is 200s. PID controllers were designed to 
compensate the system to be critically damped with a 
settling time of 15 seconds. 

Simulated temperatures are shown in Figure 2 (decay 
parameter a =■ 1 for all relevant antiwindup methods). 
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Figure 3: Closed loop temperature profiles of simulated 
furnace with continuous time anti-windup PID feedback 
laws 


AWBT method is due to the competition between the 
error e(t) and the feedback of the saturation error sig- 
nal in the AWBT signal. This can be effectively elimi- 
nated by increasing a to 100, at which point AWBT and 
VSPID are nearly indistinguishable. This behavior of the 
AWBT method is intuitively expected from the definition 
of the AWBT method; a high gain feedback of the satu- 
ration error renders the error signal e(t) inconsequential 
at the integrator summing junction. The VSPID method 
does not suffer from this drawback since the integrator is 
switched, not summed, as a function of control saturation; 
further, since the VSPID integrator settling time in satu- 
ration (approx 4 sec for a = 1) is significantly faster than 
the designed system closed loop settling time (15 sec), the 
performance of the VSPID controller is not significantly 
changed by increasing the decay factor a. 

5 Stability analysis 
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Figure 4: Closed loop temperature profiles of simulated 
furnace with continuous time anti-windup PID feedback 
laws 


We analyze the stability of the VSPID method in terms 
of a larger class of variable structure antiwindup feedback 
laws. 

Theorem 5.1 Consider a linear , time-invariant system 
S described by 

x = Ax 4- Bu + (5-2) 

y(t) = Cx(t) (5.3) 

with A = diag(A, 0), B^ ^ [ 0 0 1 ] T , = 0 

where rj is a real scalar ( integrator ). Let 

the state space be denoted as X = FR n and impose input 
saturation limits Umin> u max on the input u. Let r(t) be 
a scalar reference signal. For each K € X *, the dual 
space of X [19], define the nominal (linear) state feedback 

u n (K,x(t)) = Kx(t) and the corresponding limited state 
feedback 


and x = 


The curve labelled tf corresponds to a standard PID con- 
troller with no antiwindup law implemented. Notice that, 
for this example system, any antiwindup law is superior to 
an uncompensated PID. However, several plots are clearly 
superior to others. These signals are shown in Figures 3 
and 4. Notice that the VSPID method has reduced 
overshoot and settling time relative to the other methods 
(although CI-III sometimes appears to be competitive). 

Notice also that the VSPID and AWBT methods have 
significantly differing performance due to the nonlinear- 
ity of the VSPID method. The poor performance of the 


u s {K,x(t)) = sat(u n (if, x(£)) 5 Unxin ? Umax ) • 

For each pair (K,a) in X * x IR define the variable struc- 
ture antiwindup feedback law (VSAFL) 


u 

Ur I 


(A» 


r(t) — Cx(t) 


U — XL 5 


(5.4) 


U n 7 ^ U s 


-a(u n - w s ) 

(K, a, and x are omitted when clear by context.) Define 


X+ = {x £ X : Kx > Umax) 




X- = {x e X : Kx < Ujnin } 

X t = {x € X :x $X+L)X-} 

We say that the plant S (5.2) is in linear operation when 
x{t) £ Xt and that S is in saturation if X £ A_ U X+. 
LetV be the cone of positive definite , symmetric matrices. 
For each K £ X* define 

V(K) = {P£V : A c t P + PA C < 0, 

A c = (A + BK - B n C)} 

V{K) is the set of positive definite matrices correspond- 
ing to quadratic Lyapunov functions V(x) = x T Px for 
the system S in linear operation. For each (K,P,a,u) £ 
(X* x V x IR x LR) define the set 

V (K,P,a,u) = {x £ X : 

x T [(A - aB n K) T P + P(A - aB^K)} x 

+2 x t P(B + aB v )u < 0} 

Finally, for each (Q, 7 ) £ V x LRX , let 

X s {Qn) = {x £ X : x t Qx < 7 } C X 

be the local region in X in which we wish to stablize X. 
Then the VSAFL (5.4) with parameters K and a stabi- 
lizes the plant on X s (Q,~f) in the sense of Lyapunov if the 
following three conditions hold. 

1. P£P(K) 

2. A_ n A s C V(K, P ; a, u min ) 

3. X+r\X s CV{K,P,a : u mcix ) 

The controller globally stabilizes the plant if X s (Q,^f) = 
X S (Q, 00 ) = X . 

Proof: Observe that the sets X+, A_, Xt , V, V{K), 
V(K, P,a,u), and X 8 (P, 7 ) are all convex. The bound- 
aries of A+, A_, Xt are hyperplanes normal to the vector 
K (see Figure 5). When S is in linear operation the closed 
loop dynamics are 

x = Ax + BKx + B r} (R — Cx) 

= (.4 + BK- B n C)x + B v r (5.5) 

When the system 5 is in saturation, the closed loop dy- 
namics become 

x(t) — Ax(t) -f Bu s (t) + aB^(u^) - Kx(t)) 

= (A - aB n K) x(t) + {B + aB v )u s (t) (5.6) 



Figure 5: Sets A+, X- , and Xt for VSAFL analysis 

For each P £ V , consider the quadratic function 
V{x) = x T Px P £ V. 

A nominal state feedback u n (K ) x) stabilizes 5 if the set 
V{K) is not empty. When 5 is in saturation 

4-V(x) = X T P [(.4 - aB„K)x + (B + aB n )u s ] 

at 

+[(j4 - aB n K)x + (B + aB r) )u s ] T Px 
= x T [(>1 - aBr,K) T P + P(A - aB„K )] x 
+2 x T P{B + qB,)u s 

with u s — u m in, u max for x £ A r _ , A"+, respectively. 

Now suppose that the conditions (l)-(3) hold for a given 
set X S {Q, 7 ). Then, for every x £ X s we have that V (x) > 
0 and V{x) <0 with equality holding only at the origin. 
□ 

6 Conclusions 

We have presented the VSPID technique for the preven- 
tion of integrator windup in PID feedback control. Stabil- 
ity conditions are presented in terms of the larger class of 
VSAFL systems with state-feedback and integrator con- 
trol in Theorem 5 . 1 . All of the sets in the theorem are 
convex; further, the problem of computing a stabilizing 
nominal state feedback matrix if is a convex programming 
problem. Hence, we believe that VSAFL design problem 
can be posed as a convex programming problem[19], which 
can be solved in polynomial time [20]. However, since the 
class of VSPID controllers are not a a subset of VSAFL 
controllers, the resulting stability conditions do not lead 
to a convex programming problem. 

The combinatorial complexity of the VSAFL design 
problem increases exponentially with the number of in- 
puts subject to saturation; there are effectively 2 ^ + 1 



input laws acting in tandem, two for each input channel 
(saturation limit) and one for the linear region of opera- 
tion. 
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